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1 .  INTRODUCTION 


This  final  report  documents  work  performed  by  S&D  Dynamics,  Inc.  under 
Contract  No.  DAAK1 1 -83-C-0054  to  the  U.S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD. 

The  primary  objective  of  this  effort  is  the  extension  of  the  deterministic 
gun  dynamics  simulation  code,  DYNACODE-G,  to  accommodate  shot- to-shot  variations 
in  pertinent  forcing  functions  which  arise  during  firing,  in  order  to  assess  gun 
system  performance  on  a  probabilistic  basis.  Additional  objectives  include 
introduction  of  gun  mount  excitation  to  simulate  vehicle  motion  and  introduction 
of  an  externally  supplied  firing  signal  to  prescribe  burst-mode  fire. 

For  the  purposes  of  this  study  the  75mm  ADMAG  served  as  baseline  gun  system, 
75mm  LAV  firings  conducted  at  Aberdeen  Proving  Ground,  November  1984,  served  as 
forcing  function  data  base,  and  Ml  vehicle-hull  accelerations  for  the  Munson 
straight,  full  and  washboard  courses  served  as  vehicle  motion  data  base. 

In  order  to  achieve  the  cited  objectives,  DYNACODE-G  was  first  modified  to 
accommodate  independent  right  and  left  recoil  accumulator  loads.  Forcing  function 
data  consisting  of  propellant  generated  gas  pressure,  right  and  lvft  recoil 
accumulator  pressures  and  recoil  motion  obtained  from  seventeen  LAV  firings  were 
statistically  analyzed  by  Mrs.  Melinda  B.  Krummerich,  Mechanics  and  Structures 
Branch,  Interior  Ballistics  Division,  U.S.  Army  Ballistic  Research  Laboratory, 
to  serve  as  data  base.  DYNACODE-G  was  then  modified  to  include  a  Monte-Carlo 
routine  (with  BRL  supplied  random  number  generator)  to  operate  on  the  data  base 
to  produce  a  statistically  meaningful  sample  of  ’'mathematical"  firings.  Ml  vehicle- 
hull  accelerations  were  statistically  analyzed  by  Miss  Susan  A.  Coates,  Mechanics 
and  Structures  Branch,  IBD,  BRL,  and  DYNACODE-G  was  modified  to  accommodate  mount 
excitation  (i.e.,  simulated  vehicle  motion).  In  order  to  efficiently  handle  gun 
system  response  to  mount  excitation,  firing  excitation  and  the  response  time 
between  successive  firings,  DYNACODE-G  was  modified  by  replacing  the  fixed  time- 
step  Runge-Kutta  integration  scheme  with  an  automated  (self-determining)  variable 
time-step  routine.  Finally,  DYNACODE-G  was  modified  to  accommodate  an  externally 
sueplied  firing  signal  to  prescribe  the  timing  between  successive  rounds  in  burst¬ 
mode  fire,  based  on  gun  tube  muzzle  motion  parameters  falling  within  a  prescribed 
"window"  of  acceptable  values. 


2.  MONTE-CARLO  ROUTINE 


DYNACODE-G  was  modified  to  include  a  Monte-Carlo  routine  for  the  purpose 
of  generating  a  statistically  meaningful  sample  of  •'mathematical"  firings  based 
on  a  limited  sample  of  actual  firings.  An  overview  of  the  technique  is  depicted 
in  Fig.  1 . 

Referring  to  Fig.  1 ,  test  results  for  pertinent  forcing  function  data 
obtained  from  "M"  firings  are  statistically  analyzed.  A  random  number  generator 
is  introduced  to  produce  pertinent  forcing  function  data  for  "N"  (»  "M") 
"mathematical"  firings.  Each  set  of  forcing  function  data  obtained  for  each 
"mathematical''  firing  is  input  into  DYNACODE-G  for  the  purpose  of  generating 
muzzle  motion  data.  The  resulting  muzzle  motion  data  are  then  statistically 
analyzed  to  produce  a  probabilistic  assessment  of  muzzle  motion  for  the  class 
of  ammunition  considered. 

2 . 1  Statistical  Analysis  of  Pertinent  Forcing  Function  Data 

Forcing  function  data  consisting  of  propellant  generated  gas  pressure, 
right  and  left  recoil  accumulator  pressures  and  recoil  velocity  for  each  of 
seventeen  75mm  LAV  firings  conducted  at  Aberdeen  Proving  Ground,  November  198A, 
were  statistically  analyzed  by  Mrs.  Melinda  B.  Krummerich,  Mechanics  and 
Structures  Branch,  Interior  Ballistics  Division,  U.S.  Army  Ballistic  Research 
Laboratory. 

Individual  shot  data  for  propellant  generated  gas  pressure  was  first 
normalized  with  respect  to  time-to-peak-pressures,  t*,  as  depicted  schematically 
in  Fig.  2.  Assuming  a  normal  distribution  for  t*,  the  statistical  distribution 
for  propellant  generated  gas  pressure  as  a  function  of  time,  namely  pr(t/t*),  was 

obtained  for  the  sample  of  seventeen  firings,  as  schematically  depicted  in  Fig.  2. 

Examination  of  individual  shot  data  revealed  two  distinct  patterns  for 
right  recoil  accumulator  pressure  data,  a  single  pattern  for  left  recoil  accumulator 
pressure  data  and  a  single  pattern  for  recoil  velocity  data,  as  schematically 
depicted  in  Fig.  3.  Right  recoil  accumulator  pressure  data  for  eight  of  the 
firings,  Pattern  "A",  were  identical  with  corresponding  left  accumulator 
data;  whereas,  right  recoil  accumulator  data  for  the  remaining  nine  firings, 

Pattern  "B",  were  consistently  higher  than  corresponding  left  recoil  pressure  data. 
Statistical  distributions  for  right  recoil  accumulator  pressure  Patterns  "A"  and 
"B",  left  recoil  accumulator  pressure  and  recoil  velocity  were  obtained  as 
schematically  depicted  in  Fig.  3. 

Having  defined  the  statistical  distributions  associated  with  the  pertinent 
forcing  functions  for  the  limited  sample  of  firings,  the  Monte-Carlo  routine  was 
introduced  to  define  the  variables  associated  with  each  of  "N"  mathematical 
firings,  as  depicted  in  Fig.  A.  Referring  to  Fig.  A,  a  two-state  random  number 
generator  selects  either  Pattern  "A"  or  "B"  for  the  right  recoil  accumulator, 
while  a  normally  distributed  random  number  generator  defines  a  compatible  set  of 
forcing  function  data  for  each  of  the  "N"  mathematical  firings  for  input  (on  a 
deterministic  basis)  into  DYNACODE-G. 


Overview  of  Technique 


Fig,  3  -  Statistical  Distributions  of  Recoil  Accumulator  Pressure  and  Recoil  Motion 


Fig.  k  -  Application  of  Monte-Carlo  Routine 
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2.2  Application  of  Monte-Carlo  Routine 

In  general,  a  sample  population  size  of  thirty  data  sets  is  sufficient 
to  achieve  a  meaningful  statistical  analysis.  However,  since  a  two-state  random 
number  generator  is  needed  to  define  the  right  recoil  accumulator  pressure 
pattern,  twice  this  number,  or  sixty  data  sets,  are  appropriate.  The  Monte-Carlo 
routine  and  modified  version  of  DYNACODE-G  were  employed  as  depicted  in  Fig.  A  to 
produce  sixty,  single-shot,  "mathematical”  firings,  based  on  the  seventeen  75mm 
LAV  firings,  using  the  BRL  Cyber  7600  system.  The  results  of  these  runs  have 
been  statistically  analyzed.  It  is  noted  that  the  fixed  time-step,  fourth-order, 
Runge-Kutta  integration  scheme  has  been  retained  in  this  version  of  DYNACODE-G. 

Population  means  and  standard  deviations  were  computed  for  shot-exit  time, 
muzzle  velocity,  (total)  transverse  muzzle  displacement  at  shot-exit,  (total) 
transverse  linear  velocity  at  shot-exit,  and  (total)  pitch  rate  at  shot-exit,  as 
presented  in  Table  I.  As  may  be  seen  from  Table  I,  the  greatest  sensitivity  in 
terms  of  shot-to-shot  variations  occur  in  transverse  linear  velocity  and  pitch 
rate.  Since  these  parameters  dominate  the  "jump"  equations  it  may  be  concluded 
that  shot- to-shot  variations  can  account  for  a  10%  to  20%  increase  in  dispersion 
at  the  target. 

In  addition,  the  total  population  of  sixty  "mathematical"  firings  was 
segregated  into  two  separate  populations  characterized  by  the  distinction  between 
right  recoil  accumulator  pressure  Patterns  "A"  and  "B" .  Twenty-six  data  sets 
contained  Pattern  "A";  thirty-four  data  sets  contained  Pattern  "B".  Student  t- 
tests  on  the  resulting  means  and  standard  deviations  obtained  for  all  parameters 
examined  within  each  population  showed  no  statistically  significant  difference 
between  the  two  populations.  Hence,  it  may  be  concluded  that  the  noted  anomalous 
recoil  accumulator  pressure  patterns  have  insignificant  effect  on  the  total 
dispersion  budget. 


Finally,  a  Chi-Square  analysis  was  performed  to  determine  the  correlation 
between  transverse  linear  velocity  and  pitch  rate  at  shot-exit.  Examination  of 
the  individual  "mathematical"  firings  relative  to  the  overall  population  means 

2 

yielded  the  2x2  matrix  in  Table  II.  The  calculated  value  of  x  for  this  table, 
namely  35. A,  is  highly  significant  and  indicates  a  high  degree  of  coup]ing 
between  linear  and  angular  muzzle  velocities  at  shot-exit.  Hence,  it  may  be 
expected  that  these  parameters  will  have  an  additive  effect  in  terms  of  the 
total  error  budget  for  target  dispersion. 
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Generated  using  Monte-Carlo  technique  and  modified  version  of  DY! 
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3.  VEHICLE  MOTION 


Algorithms  were  developed  to  recover  pitch,  roll  and  vertical  accelerations 
at  t^e  gun  trunnions  in  terms  of  h>’ll  mounted  accelerometer  data  obtained  for  Ml 
vehicle  test  runs  on  the  Munson  straight,  full  and  washboard  courses.  Miss  Susan 
A.  Coates,  Mechanics  and  Structures  Branch,  Interior  Ballistics  Division,  U.S. 

Army  Ballistic  Research  Laboratory,  implemented  the  algorithms  and,  in  addition, 
performed  a  power  spectral  density  analysis  of  the  acceleration  traces. 

Sample  traces  of  pitch,  roll  and  vertical  accelerations  as  recovered  by 
the  algorithms  for  the  Munson  full  course  (BRL  Ident  1)  are  presented  in  Figs.  5a, 
6a  and  7a,  respectively.  Results  of  the  correspor ding  power  spectral  density 
analyses  are  presented  in  Figs.  5b,  6b  and  7b. 

4  functional  representation  (based  on  frequency  content)  was  developed 
to  simulate  a  typical  cycle  of  acceleration  data.  Simulated  traces  of  pitch, 
roll  and  vertical  accelerations  (based  on  the  functional  representation  developed 
and  power  spectral  density  analyses  performed)  are  presented  in  Figs.  5c,  6c  and 
7c,  respectively. 

In  addition,  expressions  for  inertia  loads  due  to  mount  excitation  were 
developed  and  incorporated  within  DYNACODE-G.  Within  the  framework  of  small 
displacement  theory,  and  in  accordance  with  the  notation  of  DYNACODE-G,  these 
loads,  applied  to  each  mass  point  of  the  gun  system,  take  the  form 
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Fig.  5a  -  IOENT  1.  bitch  acceleration 
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Fig.  6c  -  Simulated  Roll  Acceleration 
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DYNACODE-G,  with  fixed  time-step  integration  and  without  firing  excitation, 
was  exercised  first  for  the  purpose  of  demonstrating  its  ability  to  handle  gun 
system  response  to  a  typical  cycle  of  mount  excitation  (simulated  vehicle  motion) 
for  each  vehicle  course  considered.  Sample  output  obtained  for  the  Munson  full 
course  for  vertical  muzzle  displacement,  pitch,  muzzle  velocity  and  pitch  rate 
are  presented  in  Figs.  8  thru  11,  respectively.  Output  for  lateral  muzzle 
displacement,  yaw,  muzzle  velocity  and  yaw  rate  are  presented  in  Figs.  12  thru 
15,  respectively. 

Using  the  above  version  of  DYNACODE-G,  numerous  runs  were  performed  next 
to  serve  as  baseline  for  establishing  solution  convergence  of  pertinent  parameters 
as  a  function  of  integration  step-size;  a  necessary  investigation  in  the  develop¬ 
ment  of  a  variable  time-step  integration  scheme,  as  will  be  discussed  in  the  sequel. 

Although  not  a  pertinent  part  of  the  present  study,  it  is  noted  that  a 
Monte-Caric  routine,  with  equally  weighted  random  number  generator,  could  readily 
be  introduced  in  this  version  of  DYNACODE-G  should  the  user  desire  to  study  the 
stochastic  nature  of  gun  system  motion  (due  to  vehicle  motion!  just  prior  to 
firing. 
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A.  VARIABLE  TIME-STEP  INTEGRATION 


The  fixed  time-step,  fourth  order,  Runge-Kutta  integration  scheme  employed 
in  DYNACOuE-G  is  computationally  efficient  when  treating  gun  system  response  to  a 
single  or  simultaneously  applied  set  of  excitations.  In  genera] ,  the  size  of  the 
time-step  is  keyed  to  the  frequency  content  of  the  excitation  —  the  higher  the 
frequency  contei  t,  the  smaller  the  required  time-step  for  a  convergent  solution. 
However,  when  treating  gun  system  response  to  a  series  of  sequential  excitations 
with  intermittent  time  delays,  the  fixed  time-step  scheme  is  computationally 
inefficient,  particularly  in  the  time-delay  regions  where  a  far  more  coarse  time- 
step  would  prove  adequate.  To  handle  such  situations,  whose  scenario  might  consist, 
for  example,  of  tracking  gun  tube  muzzle  motion  during  firing-on-the-move,  followed 
by  a  prescribed  delay  before  the  next  firing,  computational  efficiency  dictates 
development  of  a  variable  time-step  integration  scheme. 


Using  the  solution  provided  by  PYNACOPE-G  with  fixed  time-step  integration 
for  gun  system  response  to  vehicle  motion  as  baseline,  it  was  found  that  velocity 
parameters  are  most  sensitive  to  step-size  variations,  and,  in  particular,  it  was 
found  that  convergence  of  muzzle  pitch-rate  guarantees  convergence  of  all  other 
parameters  at  all  other  mass  points;  thereby  providing  an  absolute  error  bound 
on  the  solution.  Having  identified  a  suitable  parameter  t-o  establish  convergence, 
namely  muzzle  pitch-rate,  it  remains  to  quantify  an  error  estimate  in  choosing  the 
step-size  at  a  given  time. 


The  error  produced  in  the  numerical  solution  of  a  first  order,  ordinary, 
differential  equation  using  a  fourth  order  Runge-Kutta  integration  scheme  depends 
on  the  integration  step  size.  Letting  y2  denote  the  solution  obtained  at  time 

t*  using  an  integration  step  size  of  2  At,  and  letting  y^  denote  the  corresponding 

solution  at  t*  obtained  from  two  successive  integrations  of  At  each,  the  respective 
solution  errors  are 


1 6  i  I 

“2  =  1 5  iy2  “  y1  ' 


and  e,  =-^e2. 


e2  and  c 1  serve  as  upper  and  lower  error  bounds  which  may  be  used  in  an  automated 
process  to  determine  when  to  increase  or  decrease  the  step  size,  At. 


Since  the  computation  of  requires  calculation  of  both  y1  and  y2,  an 

improved  numerical  solution  is  obtained  as  a  by-product  of  this  computation  through 
the  use  of  the  Richardson  estimator  at  tc,  namely, 


?5  71  -  h' 


It  is  noted  that  the  error  in  the  Richar.  :*on  estimator  is  a  full  order  of  magnitude 
smaller  than  the  corresponding  er^or  in  y  . 


*s?  -x-  'ns*  xzrjxr-z  j-  *a  .ysx^ar  zss  xs’a  ^rt^jrvArsxrx^JZ*  io*  xa.  r. v_ 


Based  on  the  above  considerations,  the  fixed  time-step  integration  scheme 
within  DYNACODE-G  has  been  replaced  with  an  automated  variable  tlmf -step  routine 
which  employs  the  following  logic  at  each  integration  step: 

1 .  Starting  at  time  t,  two  successive  integrations  of  step  size 
At  each  are  performed  on  the  generalized  coordinates  and 
velocities  to  determine  their  values  at  time  t  +  2  At. 

2.  Starting  again  at  time  t,  but  with  a  step  size  of  2  At, 

the  generalized  coordinates  and  velocities  are  again  determined 
at  t  +  2  At. 

3.  The  Richardson  estimators  for  generalized  coordinates  and 
velocities  are  formed  at  t  +  2  At  using  the  values  computed 
in  1  and  2.  These  estimators  serve  as  solution  at  t  +  2  At, 
as  well  as  initial  conditions  for  the  next  integration. 

A.  At  is  determined  for  the  next  integration  by  forming  upper  and 
lower  error  bounds  on  muzzle  pitch-rate  and  comparing  these 
values  against  an  acceptable  input  error,  e,  as  follows: 

If  e1  >  e,  At  is  replaced  by  At/2 
If  e2  <  e,  At  is  replaced  by  2  At 
If  e1  <  e  <  e2,  At  remains  unchanged. 

Although  the  variable  time-step  scheme  described  above  requires  one  and 
a  half  times  as  many  calculations  at  each  step  as  does  the  fixed  time-step  scheme, 
comparison  runs  show  that  convergence  to  machine  accuracy  is  achieved  with  an 
average  integration  step  which  is  sixteen  times  larger  than  the  fixed  time-step. 
Relaxation  of  the  acceptable  error,  e,  would  of  course  allow  further  improvement 
in  run  time. 
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5.  EXTERNAL  FIRING  SIGNAL 


In  order  to  optimize  the  effectiveness  of  burst-mode  fire,  the  timing 
between  rounds  of  a  burst  is  keyed  to  pertinent  muzzle  motion  parameters 
falling  within  an  acceptable  ‘'window”  of  values.  Specifically,  DYNACODE-G 
with  mount  excitation  to  accommodate  firing-on-the-move  (as  described  in 
Section  3)  and  variable  time-step  integration  (as  described  in  Section  A) 
was  modified  to  track  muzzle  displacement  and  linear  and  angular  velocities 
subsequent  to  shot-exit  of  each  round  of  a  burst.  At  each  instant  these 
quantities  are  compared  with  acceptable  predetermined  "window”  values.  An 
external  firing  signal  is  then  coupled,  at  the  user's  option,  with  either  of 
two  control  modes  —  a  vertical-plane  "window"  of  values  (without  regard  to 
horizontal-plane  motion)  and  combined  vertical  and  horizontal-plane  "windows" 
which  must  be  satisfied  simultaneously.  A  schematic  of  program  events  with 
external  firing  signal  control  is  depicted  in  Fig.  16. 

Referring  to  Fig.  16,  the  program  begins  by  initiating  vehicle  motion  at 

t=0  and  tracking  muzzie  response  in  Region  1.  At  a  preset  (but  variable)  time, 

t  ,  in  the  vehicle  motion  cycle,  a  comparison  of  muzzle  motion  parameters  versus 
s 

acceptable  "window"  values  is  initiated.  This  comparison  continues,  defining 
Region  2,  until  muzzle  motion  satisfies  "window"  values.  The  program  then 
accepts  a  firing  signal  at  tf  and  tracks  subsequent  muzzle  motion  during  the 

recoil/counter-recoil  cycle  in  Region  3.  The  program  continues  to  track  muzzle 
motion  in  Region  A,  which  accommodates  the  prescribed  time  required  to  physically 
load  the  next  round.  The  combined  time  elapsed  in  Regions  3  and  A  defines  the 
minimum  time  required  to  fire  the  next  round,  without  regard  to  satisfying 
firing  "windows".  Regions  2,  3  and  A  are  then  sequentially  repeated,  depending 
on  the  number  of  rounds  ir.  the  burst.  The  program  stops  when  the  final  round 
of  the  burst  exits  the  muzzle. 

It  is  noted  that  varying  tg  allows  the  user  to  study  the  effect  of  firing 

at  different  times  in  the  vehicle  motion  cycle.  As  previously  noted  in  Section 
3,  a  Monte-Carlo  routine,  with  equally  weighted  random  number  generator,  could 
be  introduced  to  specify  tg  on  a  probabilistic  basis. 

DYNACODE-G,  with  simulated  vehicle  motion  for  the  Munson  full  course  and 
as  herein  modified,  was  exercised  for  the  purpose  of  comparing  three  three-round 
bursts  fired  from  the  75mm  ADMAG,  with  and  without  firing  control  signals.  For 
the  purposes  of  this  illustrative  comparison,  each  three-round  burst  was 
arbitrarily  initiated  with  t  ,  and  consequently  t_,  set  equal  to  zero.  Vertical- 

S  X 

plane  "window"  values  were  arbitrarily  set  at  jy|  £  .005  in,  |vy|  £  .05  fps  and 

|w  |  -  .05  rad/sec.  Horizontal-plane  "window"  values  were  arbitrarily  set  at 

I z i  S  .01  in,  |vz|  S  .01  fps  and  |w  |  S  .01  rad/sec.  In  addition,  load  time 

was  arbitrarily  set  such  that  the  minimum  combined  time  to  cycle  and  load  is 
1  sec.  The  user  may  of  course  vary  t  ,  "window"  values  and  load  time  according 
to  specific  needs.  ° 


Time-histories  of  vertical  muzzle  displacement,  velocity  and  pitch  rate 
are  presented  graphically  in  Figs.  17  thru  19  for  a  three-round  burst  without 
firing  control  signal,  in  Figs.  20  thru  22  with  a  vertical-plane  "window"  and 
in  Figs.  23  thru  25  with  combined  vertical  and  horizontal-plane  "windows". 
Corresponding  time-histories  of  lateral  muzzle  displacement,  velocity  and  yaw 
rate  are  presented  graphically  in  Figs.  26  thru  28  for  the  three-round  burst 
without  firing  control  signal,  in  Figs.  2y  thru  31  with  a  vertical-plane  "window" 
and  in  Figs.  32  thru  34  with  combined  vertical  and  horizontal-plane  "windows". 

In  addition,  a  summary  of  pertinent  muzzle  motion  parameters  at  shot-exit  of 
each  round  of  each  burst  is  presented  for  comparison  purposes  in  Table  III . 

As  may  be  seen  from  the  figures  and  comparisons  presented,  it  is  not 
readily  obvious  that  requiring  individual  muzzle  motion  parameters  prior 
to  firing  be  bounded  within  acceptable  "window"  values  guarantees  optimal 
muzzle  motion  at  shot-exit,  particularly  when  firing  on  a  moving  vehicle. 

The  concept  of  a  firing  "window"  indeed  seems  plausible;  however,  it  is  not 
yet  known  which  parameters  or  combinations  require  bounding.  In  addition, 
it  is  rather  obvious  that  when  firing-on-the-wiove,  consideration  must  also  be 
given  to  the  nature  of  vehicle  motion,  the  relative  timing  of  the  firing 
signal  and  phase  relations  between  muzzle  motion  and  vehicle  motion  prior 
to  firing. 


Vertical  Muzzle  Displacement  for  Three-Round  Burst  Without  Firing  Control 
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Fig.  17  -  (Cont'd 


Fig.  18  -  Vertical  Muzzle  Velocity  for  Three-Round  Burst  Without  Firing  Control  Signal 
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Fig.  19  (Cont’d. 
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Fig.  20  -  Vertical  Muzzle  Displacement  for  Three-Hound  Burst  With  Vertical-Plane  '’Window 
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Fig.  20  (Cont'd. 


Vertical  Muzzle  Velocity  for  Three-Round  Burst  With  Vertical-Plane  "Window" 
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Fig.  2k  -  Vertical  Muzzle  Velocity  for  Three-Round  Burst  With 
Combined  Vertical  and  Horizontal- Plane  "Windows" 
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Fig.  25  -  Muzsle  Pitch  Rate  for  Three-Round  Burst  With 

Combined  Vertical  and  Horizontal-Plane  "Windows” 
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Fig.  26  -  Lateral  Muzzle  Displacement  for  Three-Round  Burst  Without  Firing  Control  Signal 
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Fig.  27  (Cont'd. 
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Fig.  29  -  Lateral  Muzzle  Displacement  for  Three-Round  Burst  With  Vertical-Plane  "Window*1 
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Fig.  32  (Cont'd. 
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Fig.  33  (Cont’d 
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Without  firing  control  signal 

Vertical-plane  "window”  only 

Combined  vertical  and  horizontal-plane  "windows 


6.  CONCLUSIONS 


Employing  a  Monte-Carlo  routine  and  random  number  generator,  DYNACODE-G 
has  been  modified  such  that  based  on  a  limited  number  of  actual  firings  the 
user  may  now  generate  a  statistically  meaningful  sample  of  "mathematical" 
firings  in  order  to  assess  shot-to-shot  ar '.ations  in  muzzle  motion  (or  any 
other  gun  system  response  parameter  of  ,t  rest),  for  a  particular  class^pf 
ammunition.  Although  the  75mm  ADMAG  .«■*•  >■"  d  as  baseline  gun  system,  with 
LAV  firings  as  data  base,  the  techniqi'  developed  are  general  and  may  be 
applied  to  any  other  gun  system  and  data  base  of  interest. 

To  simulate  gun  system  response  to  vehicle  motion  over  prescribed  terrain, 
DYNACODE-G  has  been  modified  tc  accommodate  acceleration  inputs  at  the  trunnion 
level.  Using  the  75mm  ADMAG  as  baseline  gun  system,  and  Ml  vehicle  motion  on 
the  Munson  straight,  full  and  washboard  courses  as  data  base,  typical  cycles 
of  vehicle  motion  data  were  simulated  analytically  and  input  into  DYNACODE-G. 
Resulting  output  provides  a  means  for  assessing  the  effects  of  firing-on-the- 
move  and  the  user  may  readily  introduce  a  Monte-Carlo  routine,  with  equally 
weighted  random  number  generator ,  to  predict  probable  gun  tube  muzzle  motion 
prior  to  firing  a  vehicle  mounted  gun.  The  analyses  performed  are  restricted 
to  the  vehicle  motion  data  base  provided;  however,  the  techniques  used  in 
simulating  vehicle  motion  should  prove  generally  applicable  to  other  vehicles 
and  terrain. 

In  order  to  efficiently  handle  the  comDutatienal  requirements  of  tracking 
gun  system  response  during  burst-mode  fire,  an  automated  (self-determining) 
variable  time-step  numerical  integration  scheme  has  been  developed  and  success¬ 
fully  implemented.  This  scheme  could  of  course  be  introduced  into  the  other 
versions  of  DYNACODE-G:  however,  the  fixed  time-step  routine  has  proven  tc  be 
of  adequate  computational  efficiency  when  handling  single-shot  f-re. 

Finally,  a  version  of  DYNACODE-G  with  variable  time-step  integration 
and  vehicle  motion  has  been  extended  to  incorporate  an  externa]  firing  signal 
and  firing  "windows"  for  controlled  burst-mode  fire.  The  program  logic  developed 
tracks  muzzle  motion  subsequent  to  a  firing  and  accept?  a  firing  signal  for  the 
next  round  of  a  burst  after  sufficient  time  has  elapsed  (i)  to  load  the  next 
round  and  (ii)  for  muzzle  motion  to  subside  to  acceptable  "window"  values. 

Using,  for  illustrative  purposes,  a  loau  time  of  1  sec  for  the  75mm  ADMAG., 
muzzle  motion  has  generally  dissipated  when  firing  on  a  fixed  mount;  however, 
muzzle  motion  is  reinforced  when  firing  on  a  moving  vehicle.  Using  simulated 
Ml  vehicle  motion  on  the  Munson  full  course,  three  three-round  bursts,  consisting 
of  uncontrolled  fire,  controlled  fire  with  a  vertical  plane  "window"  only,  and 
controlled  fire  with  both  vertical  and  horizontal  plane  "windows"  have  been 
'■ompared.  The  results  of  this  '■'cmparison  show  that  although  it  indeed  seems 
plausible  to  introduce  a  firing  "window"  to  optimize  burst-mode  fire,  further 
study  is  needed  to  determine  the  parameters  (and/or  combinations  thereof)  which 
must  be  bounded,  as  well  as  acceptable  "window"  values.  The  latter  *111  of 
course  require  a  trade-off  with  delay  time  between  firings  to  achieve  burst-mode 
optimization.  In  addition,  when  f iring-on-the-move,  consideration  must  also  be 
given  to  the  nature  of  vehicle  motion,  the  relative  timing  of  the  firing  signal 
and  phase  relations  between  muzzle  and  vehicle  motions. 
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